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ABSTRACT 

A general framework for tests of Lorentz invariance with electromagnetic waves is presented, al- 
lowing for operators of arbitrary mass dimension. Signatures of Lorentz violations include vacuum 
birefringence, vacuum dispersion, and anisotropies. Sensitive searches for violations using sources such 
as active galaxies, gamma-ray bursts, and the cosmic microwave background are discussed. Direction- 
dependent dispersion constraints are obtained on operators of dimension 6 and 8 using gamma-ray 
bursts and the blazar Markarian 501. Stringent constraints on operators of dimension 3 are found 
using 5-year data from the Wilkinson Microwave Anisotropy Probe. No evidence appears for isotropic 
Lorentz violation, while some support at la is found for anisotropic violation. 

Subject headings: relativity — gravitation — cosmic microwave background — gamma rays: bursts 
— galaxies: active 



Recent years have seen a resurgence in tests of rel- 
ativity, spurred in part by the prospect of relativ- 
ity vi olations arising in a unif i ed description of na- 
ture ([Kostelecky fe Samuel! 119891 : IKostelecky fc Potting! 
1991). Experimental searches for violations of Lorentz 
invariance, the symmetry underlying relativity, have 
been perfo rmed in a wide range of sy stems (for data 
tables, see IKostelecky fe Russelll ((2008)). Historically, 
experiments probing the behavior of light have been 
central in confirming relativity. Contemporary ver- 
sions of the classic M ichelson-Morley and Kennedy- 
Thorndike experiments (fLipa et al"1 l2b03: Ant onini et all 
12001 iMuller et al.l 12007ft remain among the most sensi- 
tive tests today. 

Some tight constraints on relativity violations 
have been achieved by seeking tiny changes in 
light that has propagated over astrophysical dis- 
tances. Many of these search for a change in 
polarization resulting from vac uum birefringence , 
using sources such a s galaxi es (ICarroll et al.l 119901: 
Colladav fc Kostelecky! Il998t IKostelecky fc Mewes! 
2001L 12002ft gamma -ray bursts (GRB) (iMitrofanoyl 



2003: Uacobson et all [200l IKostelecky fc Mewes! 1200 



Kahniashvili et alj 120061 : 



Fan et al 



cosm i c microwave back ground 
2006: iGamboa et all [200 
Cabella et al.l l2007t 



lOQj and the 
ound (CMB) (iFeng et~ap 
IKostelecky fc Mewes! 120071: 
120091 : iXia et ail 



Komatsu et al 



2008tlKahniashvili et al.ll2008f ). Others seek a frequency- 



dependent velocity arising from vacuum dispersion, us- 
ing GRB, pulsars, and blaz a rs (lAmelino-Camelia et all 
19981: IKostelecky fc Mewesl 120021: iBoggs et al.l 12004 



Martinez fc Piranl 120061: lEllis et all (20061 : lLamon et all 
20081 : [Albert et alJl2008f l. Here, we present a general the- 
oretical framework that characterizes Lorentz-violating 
effects on the vacuum propagation of electromagnetic 
waves and includes operators of all mass dimensions. We 
discuss several techniques that can be used to search for 
the unconventional signals of Lorentz violation, Using 
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vacuum-dispersion constraints from GRB and the blazar 
Markarian 501, we place new direction-dependent limits 
on several combinations of coefficients for Lorentz vio- 
lation. We also perform a search for Lorentz violations 
in the 5-year results from the Wilkinson Mic r owave 
Anisotropy Probe (WMAP) (iKomatsu et all 120091 : 
iHinshaw etaLl 120091: iNolta et all 12009! ) . finding some 
evidence for anisotropic violations but no support for 
isotropic violations. 

At attainable energies, violations of Lorentz invari- 
ance are described by a fr amework called the Standard- 
Model Extension (SME ) (jColladav fc Kostelecky! 119971 . 
Il998t IKostelecky! 12004ft that is bas ed on effective field 
theory (|Kosteleckv fc Potting! 1 1995ft. Approache s out- 
side fieldl!li«i>ry~als^^ The 
SME characterizes all realistic violations affecting known 
particles and fields, while incorporating otherwise estab- 
lished physics. Much of the work on Lorentz violation has 
focused on the minimal SME, which restricts attention 
to gauge-invariant operators of renormalizable dimen- 
sion. In this work, we consider the gauge-invariant pure- 
photon sector of the full SME with Lorentz-violating op- 
erat ors of arbitrary dimension, w hich has Lagrange den- 
sity ([Kostelecky fc Mewesl [200% 

C = -\F ilv F^ + ±e KX ^A x (kAF).F^ 

-\F KX {k F T x ^F^ , (1) 

where A M is the 4-potential with field strength F^. In 
a flat background with energy-momentum conservation, 
the Lorentz violation arises through the differential op- 
erators 



(d) , a l--- a (d-3) 



d ai . . . d a(d 3) , (2) 



d— odd 



(k F r x ^=^2 {kfr x ^- a ^d ai ...3 Q(d _ 4) .(3) 

d— even 

The constant coefficients (k^ F ) K (d 3> and 
{k^) KX ^ ai --- a {d-i) characterize the degree of Lorentz 
violation. The former control CPT-odd operators and 
are nonzero for odd dimension d > 3, while the latter 
control CPT-even operators and are restricted to even 
d> 4. 



2 



The Lagrange density |T]) yields modified Maxwell 
equations. At leading order in coefficients for Lorentz vi- 
olation, two plane- wave solutions exist. The correspond- 
ing two modified dispersion relations can be written in 
the form 



pH^i^^vW+FF+Wk (4) 

where p and w are the wavenumber and frequency, re- 
spectively. It follows that electromagnetic waves gener- 
ically contain two propagating modes with different ve- 
locities and polarizations. The symbols <; , q , <j 2 , and <j 3 
represent certain combinations of coefficients for Lorentz 
violation, and they depend on the frequency u> and direc- 
tion of propagation p. With convenient normalizations, 
s 2 , and <r 3 are the Stokes parameters s 1 = Q, s 2 = U, 
and s 3 = V of the faster mode, while is a scalar com- 
bination analogous to the intensity s = I. These four 
combinations completely control the leading-order effects 
of Lorentz violation on light propagating through empty 
space. The combination <; 3 depends only on the coeffi- 
cients (k^p) K {d 3> , while q°, s 1 , and <r 2 depend only 
on the coefficients (fc^)"*""" 1 - ^-*). 

It is convenient to identify a minimal set of coeffi- 
cient combinations that affect light propagating in vacuo. 
This can be accomplished through spherical-harmonic 
decomposition. Since <; , <; 3 are rotation scalars while 
? , <r 2 are rotation tensors, their decomposition must in- 
volve some form of tensor spherical harmonics. The spin- 
wei ghted harmonics s Yj m (p) provide a well- u nderstood 



set ((Newman k Penrosdll966t lGoldberdll967l) . The in- 
dex s is the spin weight, which up to a sign is equivalent 
to helicity. Decomposing yields 

djm 



djm 
djm 



(B)jm) ' 
(V)jm ' V>) 



where j < d — 2 and h = —p is a unit vector pointing to 
the source in astrophysics tests. 

With this decomposition, all types of Lorentz viola- 
tions for propagation in vacuo can now be simply char- 
acterized using four sets of spherical coefficients, k^ d y m , 

k \% jm > k % m fOT CPT-even effects and k$ )jm for CPT- 
odd effects. For each coefficient, the underlying Lorentz- 
violating operator has mass dimension d and eigenvalues 
of total angular momentum given by jm, as usual. For 
light from astrophysical sources, dispersion arises when 
the speed of propagation depends on frequency, which oc- 
curs for any nonzero coefficient with d ^ 4. Birefringence 
results when the usual degeneracy among polarizations is 

broken, for which at least one of kfj )jm , kf^ )jm , fc ( ( y )jm 
is nonzero. For example, all operators producing light- 
speed corrections that are linear in the energy have d = 5 
and are necessarily birefringcnt. The only coefficients for 
nonbirefringent dispersion are therefore k^- m with even 

d > 6. Since birefringence tests using polarimetry are 
typically many orders of magnitude more sensitive than 
dispersion tests using timing, in the following discussion 



of dispersion we focus only on coefficients for nonbire- 
fringent dispersion. 

Tests for vacuum dispersion seek differences in the ve- 
locity of light at different wavelengths. In the present 
context with zero bircfringent coefficients, the change 
in velocity is Sv ~ — We see from Eq. ([5]) that 
the velocity generically depends on the direction n as 
well as the frequency to. Typical analyses study ex- 
plosive or pulsed sources of radiation producing light 
over a wide wavelength range in short time periods, 
comparing the arrival times of different wavelengths. 
This idea has been the focus of many searches based 
on m o dified dispersion relati o ns (lAmelino-Camelia et al.l 
1 19981: iKosteleckv k Mewesl l2002t iBoggs et all 12004 
Martinez k Piranl l2006t lEllis et all (20061 : iLamon et al l 
l2008t I Albert et al.l 120081 ). Many of these studies as- 
sume isotropic violations, which corresponds to the limit 
j = m = 0. However, at each dimension d, this isotropic 
restriction misses (d 2 — 2d — 2) possible effects from 
anisotropic violations. 

To calculate arrival-time differ ences in an expandin g 
universe, some care is required ((Jacob k Piranl [20081 . 
In the present case, the photons propagate between two 
comoving objects, so the relevant coordinate interval is 
dl c = (1 + z)dl p = —v z dz/H z . Here, v z is the parti- 
cle velocity at redshift z, and H z = i/o(^rC 4 + ^mC 3 + 

n k ( 2 + n A y/ 2 with c = 1 + z is the Hubble expansion 

rate at z in terms of the present-day Hubble constant 
Ho ~ 71 km/s/Mpc, radiation density fl r ~ 0, matter 
density fl m ~ 0.27, vacuum density fljy — 0.73, and cur- 
vature density O/j = 1 — Q r — Q m — The total coor- 
dinate distance is the same for all wavelengths, but the 
travel times may differ. Integrating dl c from the same 
initial time to the two arrival times for the two velocities 
gives a relation for the arrival-time difference At, which 
depends on the two energies and the source location on 
the sky. For the present case with Lorentz violation at 
dimension d, we find 



At 



-Aw 



d-4 



(1 + *) 



d-4 



jm 



where Aoj d 4 is the difference in uo d 4 between the two 
frequencies. 

As an illustration, consider the bright gamma-ray burst 
GRB 021206 at right ascension 240° and declination 
—9.7°. Over energies from 3 to 17 MeV, arrival-time dif- 
ferences arenajnore than At < 4.8 ms for this source at 
z ~ 0.3 (|Boggs et al.ll2004( ). Numerical integration of Eq. 
([6]) leads to a bound on one direction-specific combina- 
tion of the 25 independent coefficients for nonbirefringent 
dispersion with d = 6: 



E 



y 3m (99.7 o , 240°) k$ jm < 1 x 10- 16 GeV- 2 . (7) 



For the 49 independent nonbirefringent dispersive oper- 
ators with d = 8, we obtain 



^ y 3m (99.7 o ,240 o )^ m < 3 x 10- 

jm 



GeV" 



(8) 



Operators with higher d can be treated similarly. Note 
that many sources are required to constrain fully the 
coefficient space for a given d. In contrast, only one 



3 



source is needed to constrain fully the corresponding co- 
efficient in the restrictive isotropic limit j = m = 0. In 
this limit, the bounds and ((SJ reduce to k?j\ 00 < 

4 x 1(T 16 GeV" 2 and fc ( ( fj 00 < 9 x 1CT 13 GeV~ 4 , respec- 
tively. 

As another example, consider Markarian 501, which 
lies at z ~ 0.03. This source produces flares with photon 
energies in the TeV range, making it particularly sensi- 
tive to an en ergy-dependent velocity and also to thresh- 
old analyses (|Amelino-Camelia fc Pira"r] [2001). A recent 
analysis of observations by the MAGIC collaboration 
found some evidence for a nonbircfringent v elocity defect 
of th e form Sv = -uj/M or Sv = -uj 2 /M 2 (|Albert et all 

2008) . The first case is incompatible with the present 
treatment; a reanalysis incorporating the necessary bire- 
fringence could yield comparatively weak but compatible 
new bounds. The second case suggests dispersion with 
M ~ x 10 10 GeV, assuming an arrival-time lag due 
entirely to nonbirefringent Lorentz violation. For d = 6, 
this yields the single constraint 

^ y jm (50.2°, 253°) fc ( % m ~ x 10~ 22 GeV" 2 , 

jm 

(9) 

consistent with the GRB bound ||7J). In the isotropic 
limit, this becomes kfl m ~ 10 x 10" 22 GeV" 2 . 

Next, we consider tests for vacuum birefringence. In 
birefringent scenarios, the two plane- wave eigenmodes 
travel at slightly different velocities, which alters their 
superposition and hence the net polarization of the light 
as it propagates in free space. The polarization change 
is equivalent to a rotation of the Stokes vector s = 
(s 1 ,s 2 ,s 3 ) T about the birefringent axis <; = (c 1 , <T 2 , <T 3 ) T . 
The total rotation angle is equal to the relative phase 
change between the two eigenmodes. Infinitesimally, 
the rate of rotation is ds/dt — — iS ■ s, where S ab = 
— 2ioje abc s c is the rotation generator. Integration from 
source redshift z to taking into account the cosmologi- 
cal expansion yields the net change in the Stokes vector, 

A3 = f ;l z ■ * dz , (io) 

where S 2 is the rotation matrix at the blue-shifted fre- 
quency (1 + z)u> and source direction n. The net polar- 
ization change As can depend on frequency and direc- 
tion of propagation. To seek birefringence, we can either 
model the polarization at the source and search for dis- 
crepancies in the observed polarization, or we can test 
for unexpected frequency dependences. 

In what follows, we investigate vacuum birefringence 
via the CMB, leaving the use of G RB polarimetry in this 
conte xt to be discussed elsewhere ()Kosteleckv fc Mewed 

2009) . The CMB has a long baseline but comparatively 
low frequency, which implies lesser sensitivities to d > 3 
violations relative to higher-frequency sources. Here, we 
focus on the four d = 3 Lorentz-violating operators. 
These induce energy-independent polarization changes, 
so the best constraints are expected from the most dis- 
tant sources irrespective of frequency. The CMB there- 
fore has the potential to yield maximal sensitivity to 
these CPT-odd operators. For any CPT-odd case, bire- 
fringence causes a rotation of the Stokes vector about 



the s 3 axis, corresponding to a rotation of the linear- 
polarization angle ip with no change in the degree of lin- 
ear or circular polarization. For d = 3, the value of tj) 
at present is ip = tp z + 5ip z , where ipz is the blueshifted 
angle and 5tp z is its rotation, 

**.=}' ^|^EoW*)*fV (id 

u \ / j m 

Taking z = 1100 for the CMB and including a small 
radiation component Q r ~ 0.015, the rotation reduces to 
the direction-dependent approximation 

^cmb 3.5° x 10 43 GeV-^oVjmlfl) k$ )jm . (12) 

jm 

We remark in passing that CPT-even operators produce 
a complicated mixing of circular a nd linear polarization, 
rathe r than a simple rotation of ip (|Kosteleckv fc Mewesl 
f2009h . 

We next search for the above effect in the recent 
WMAP 5-year results (jHinshaw et alj 120091 ; iNolta et all 
2009). We generate initial sky maps of the Stokes param- 
eters using the best fit correlation coefficients Cj as cal- 
culated by the WMAP collaboration within the A-CDM 
model assuming gaussianality. The Stokes parameters at 
each point on the sky are then rotated appropriately and 
used to calculate the Cj coefficients at the present epoch. 
The likelihood of these coefficients is determined using 
available WMAP software. The underlying cosmology is 
kept fixed, so we are comparing the likelihood of Lorentz 
violation relative to a reasonable Lorentz-invariant cos- 
mology. Our analysis uses TE and TB data at high-Z 
corresponding to j = 24-450, disregarding the TT data. 
The latter is a good approximation because the TT data 
would dominate an analysis with varying cosmology and 
therefore hold the cosmology comparatively fixed. 

The correlation coefficients Cj are rotationally invari- 
ant, so our analysis has sensitivity only to rotationally 
invariant combinations of Lorentz-violating coefficients. 
In the present context, these are the isotropic coeffi- 
cient (k^p) T = — fc^y) 00 /'\/47r and the scalar magnitude 

\k ( 2l\ = (6|^ )n | 2 + 3|^ )10 | 2 ) 1/2 /^. Inparticu- 

(3) 

lar, our results are independent of the direction of k\ F . 
Although the analysis contains no a priori anisotropies, 
the procedure involves generating random realizations 
that contain anisotropies. As a result, the likelihood 
L(fc/yw m ,r) for a given realization r is anisotropic. In 

obtaining the total likelihood for a given set of coefficients 
for Lorentz violation, we sum over the likelihoods of all 
possible realizations weighted by the probability density 

P(r), yielding L(k$ )jm ) = £ r P(r)i(ft m ,r) . This 
total likelihood is simply the average over all possible uni- 
verses and is a rotationally invariant indicator of Lorentz 

violation. Here, we estimate L(fc|y^ m ) for a range of 

values of fejU^ by averaging over 3,000 realizations per 
value. The results for the four coefficients with d = 3 are 
shown in Fig. [TJ 
For the isotropic coefficient, we obtain the la result 

k< (v)oo = ( 2 - 3 ± 5 - 4 ) x 1°" 43 GcV - (I 3 ) 
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and k 



(3) 



1,(3) 
(V)00 




„„, - (2.5 ± 3.0) x 10~ 43 GeV (|Cabella et alJ 

l2007f ) . Some indication of a nonzero rotation has been 
found in pre vious studies. One involving BOOMERANG 
(B03) data (jMontrov et al.|[200l [pTacentini et al.ll2006l) 

alone yielded the possibilit y fc ( ( ^ )00 = (12 ± 7) x 1(T 43 

GeV (|Kosteleckv fc Mewesll2007f) . while the result from 
another analysis combining B03 and WMAP 5-year 



data corresponds to fc[^ )00 < (2.6 ± 1.9) x 1CT 43 GeV 

(|Xia et al.l 120081 ). While consistent with the latter, our 
result (|13p shows little evidence for isotropic Lorentz vi- 
olation. 

For each anisotropic coefficient, Fig.Q]displays the like- 
lihood. As expected, the results are independent of the 

direction of k AF . The plot symmetries reflect the ex- 



pected behavior under k AF 
the result 



(3) 

k AF . These plots yield 



Im k. 



(V)ll 



Fig. 1. — Relative likelihood versus the four d = 3 coefficients for 
Lorentz violation. Points indicate the values at which the ensemble 
averages were made, and the bars represent their standard errors. 
The line is an extrapolation through the points. Dark-gray regions 
indicate the estimated 68% confidence interval, while the light-gray 
region shows the 95% level. All coefficients are in units of 10~ 42 
GeV. 

This improves by about an order of magnitude on 
the p revious limit from radiogalaxy data (| Carroll et al.l 
1990). The result is also consistent with that of 
the WMAP collab oration, which found a rotation of 
Sip ~ 1.2° ± 2.2° ()Komatsu et al.l l2009h corresponding 

to fc^Ln < (1-2 ± 2.2) X 10" 43 GeV, and that of an- 



.00 
"AF 



= (10t 4 ) x 10~ 43 GeV, 



(14) 



V (V)00 



(3) 

other recent analysis yielding k)y^ 00 < 2.5 x 10 

(jKahniashvili et all l2008f ) . Othe r reported limits ob- 
tained from WMAP 3-year data (jPage et all I2007T ) in- 
clude fc^ = (6.0 ±4.0) x 10~ 43 GeV (|Feng et ali\2~(M) 



-43 



GeV 



revealing some evidence at la for anisotropic Lorentz 
violation in the WMAP 5-year data. This agrees with 
the indication of anisotropic Lorentz violation found from 

( 3) 

an analysis of B03 data, which corresponds to |fe^p| = 
(15 ± 6) x 10~ 43 GeV (IKostefeckv fc Mewesl[2?j07l) . The 
data are consistent with no Lorentz violation at 2<r, with 

a 95% confidence level of \k ( AF \ < 2 x 10~ 42 GeV. This 

(3) 

fully constrains the vector components of k AFl and the 
results (TT3"|) and (fl4"|) provide a measurement of all four 
of the d — 3 coefficients for Lorentz violation. 
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